Merging of single-particle levels in finite Fermi systems 
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Properties of the distribution of single-particle levels adjacent to the Fermi surface in finite Fermi 
systems are studied, focusing on the case in which these levels are degenerate. The interaction of 
the quasiparticles occupying these levels lifts the degeneracy and affects the distance between the 
closest levels on opposite sides of the Fermi surface, as the number of particles in the system is 
varied. In addition to the familiar scenario of level crossing, a new phenomenon is uncovered, in 
which the merging of single-particle levels results in the disappearance of well-defined single-particle 
excitations. Implications of this finding are discussed for nuclear, solid-state, and atomic systems. 
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Advanced technologies seek to exploit the properties 
of objects of nanometer size in the design of new ma- 
terials and devices. The behavior of the electronic sys- 
tem within such a nanoscale object is largely determined 
by the structure of individual single-particle (sp) levels. 
It therefore seems opportune to revive the fundamental 
study of single-particle aspects of finite Fermi systems 
as developed many years ago for atomic nuclei^ in the 
expectation that the findings may also be of value for 
electronic systems of current technological importance. 

In homogeneous matter, all the relevant measurable 
quantities, such as various susceptibilities, arc functions 
of a single momentum transfer variable q. Inhomoge- 
neous systems with a uniform distribution of sp levels 
possess basically the same properties. However, the situ- 
ation changes when the spectrum of their sp excitations 
is degenerate. This degeneracy implies the existence of a 
new energy scale -D mm , given by the difference between 
the energies of the closest sp levels lying on opposite sides 
of the Fermi surface. The properties of such systems ex- 
hibit striking departures from what is found in homoge- 
neous matter. To offer a prominent example, consider the 
ground-state energy Eq (A) of atomic nuclei as a function 
of mass number A. For most nuclei, this quantity is well 
described by the Bethe-Weiszacker liquid-drop formula. 
However, nuclei with a magic number of protons or neu- 
trons have spherical form, and the relevant energy scale 
D mm is several times larger than the average distance 
between neighboring sp levels in non-magic nuclei. This 
energy spacing provides a shell correction 5E S , lowering 
the liquid-drop binding energy and rendering the ground 
states of known magic nuclei stable with respect to any 
mode of decay. 2 Another example is associated with the 
degeneracy of the sp spectrum of the two-dimensional 
electron gas in an external magnetic field. In this case, 
the degeneracy gives rise to a step-like behavior of the 
chemical potential fJ>(A), triggering oscillations of ther- 
modynamic properties^ 

Customary explanations of such extraordinary behav- 
ior do not take account of the alteration of key quantities 
due to interactions between added particles under varia- 



tion of their number. In many cases such an approxima- 
tion is justified, since these interactions do not affect the 
deviations mentioned above, even if the sp levels cross one 
another. However, we shall demonstrate that the depen- 
dence of sp energies e s — E s (A+1)—Eq{A) on quasiparti- 
cle occupation numbers ri, inherent in Fermi-liquid (FL) 
theory, allows for an alternative scenario. The familiar 
crossing of sp levels is replaced by a merging of these lev- 
els, a new phenomenon that leads to the disappearance 
of well-defined sp excitations and drastic departures from 
predictions of standard FL theory. 

To gain insight into this unconventional scenario, we 
first consider a schematic model involving three neutron 
levels in an open shell of a spherical nucleus. The levels 
are denoted — , 0, and +, in order of increasing energy, 
and the distance between — and and between and + 
has the same value D. As usual, the sp energies and wave 
functions <p\(r) = Rni{ r )^ jimi'n) as well, are solutions 
of equation [p 2 /2M + T,(r,p)]ip\(r) — e\(p\(r), where E 
stands for the self-energy. In even-even spherical nuclei, 
which have total angular momentum J = due to pairing 
correlations, the energies e\ are independent of the mag- 
netic quantum number m associated with the total sp an- 
gular momentum j. We suppose that the level — is filled, 
the level + is empty, and N neutrons are added to the 
level 0, changing the density by Sp(r) = N R^ olo (r) / An . 
We assume that l_^l ^l + ~ A 1 / 3 > 1. 

It is our goal here to explore the consequences of the 
dependence of the sp energies t\(n) on the distribution n. 
In what follows, we shall retain only a major, spin- and 
momentum-independent part of the self-energy E and a 
primary, 5(r)-like portion of the Landau-Migdal interac- 
tion function /. Accordingly, the FL relation between E 
and the density p responsible for the variation of e\(n) 
with n is simplified toi <5E(r) = f[p(r)]Sp(r). For the 
sake of simplicity, the diagonal and nondiagonal matrix 
elements of / are assigned the respective values 

u = |i? 2 ; (r)/[p(r)]i? 2 ; (r)r 2 dr/47T, (1) 
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R 2 nl (r)f{p(r)]R 2 nill (r)r 2 dr/4ir, (2) 



independently of the quantum numbers nl, ri\l\. 

A simple estimate of the ratio u/w is obtained using 
a semiclassical approximation R n i(r) ~ r" 1 cos J" p(r)dr, 
with the result u ~ 3w/2. We next observe that /(p) 
is positive at densities close to equilibrium^ but changes 
sign as p — > 0; hence the signs of u and w may depend 
on the quantum numbers of the sp levels in play. 

Based on these results, the dimensionless shift Sk(N) — 
[e k (N) - e fc (0)] /D for fc = 0, +, - is given by 

e (N)=n U 1 e+(N) = e_(JV) = n W , (3) 

where = Nk/(2jk + 1) is the occupation number of 
level k, U = u(2j + 1)/D, and W = w{2j + 1)/D. It 
is readily verified that if fppM/-K 2 ~ 1, where pp — 
\f2Mep and ep is the Fermi energy, then the integral 
(0) has a value u ~ £f/A and therefore |J7| ~ 1, since 
D ^ ep /A 2 / 3 in spherical nuclei, 

According to Eqs. J3J, the distance e+(N) — e-(-W) 
remains invariant when iV increases. On the other hand, 
the difference d+(N) = £+(N) — eo(N) decreases with 
N when U > W > 0, as does the distance d-(N) — 
e (N) — e-(N) in the opposite case, U < W < 0. 

Now let us determine what can happen when the func- 
tions d±(N) change their signs before the sp level is 
completely filled. We first examine the case U < W < 0. 
According to Eqs. (0, the sign of d-(N) changes at 
ri 0c = 1/(W — U), which requires W — U to be greater 
than 1 to meet the restriction n < 1. The usual Hartree- 
Fock (HF) scenario prescribes that for uq > uq c , quasi- 
particles must leave the occupied level — and resettle 
into the unfilled level 0. Further, when the dependence 
e\ = e\(n) from Eqs. J3J is brought into the picture, this 
effect is seen to promote the HF rearrangement. 

In the opposite case, U > W > 0, the function d+(N) 
changes sign at n 0c = 1/(U — W), implying U — W > 1. 
In order to satisfy this inequality, the repulsive part of 
the interaction / has to be sufficiently large, or else the 
scale D must be rather small. At uq > tiq c , the HF 
scenario requires the quasiparticlcs to leave the unfilled 
level and move into the empty level +. Were this sce- 
nario the correct one, the rearranged sp energies would 
obey the equations €q{N) — eo(N c ) + SN c (w — u) and 
e+(N) = eo(N c )+SN c (u—w), where SN C is the number of 
quasiparticles shifted from level to level +. The S term 
in each of these equations arises due to the feedback of the 
immigrating quasiparticles. Upon subtracting one equa- 
tion from the other, we find that e+(N) — €q(N) > for 
any SN C > 0, which says that the level + lies above rather 
than below the level 0. We thus arrive at a contradiction 
that excludes the HF scenario in the case U > W > 0. 

Under these conditions, a new ground state must form, 
denoted henceforth by M. As will now be shown, in the 
state M both of the levels and + are partially occupied. 
Solution of the problem for this case reduces to finding 



the minimum of the relevant energy functional 



E = e Q (0)N Q + e+(0)N+ + - [u{N 2 + Nf.) + 2wN N+\ 

. (4) 

with Nk — X)m n km, through the variational conditions 
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where p is the chemical potential. Such a condition for 
characterization of a rearranged ground state first ap- 
peared in Refs. 0, where homogeneous Fermi systems 
were addressed, without attention to degeneracy of sp 
levels. Eqs. (JHJ) are conveniently rewritten as a condi- 
tions for the coincidence of the sp energies eg and e + , 



eo(JV) = eo(0) + N u + N+w = /i , 
e+(N) = e+(0) + N w + N + u = fi , 

which, at N > N c = (2j + 1)/(Z7 - W), yield 
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FIG. 1: Top panels: Dimensionless distance d+ = (e+—eo)/D 
between levels + and as a function of the ratio x — N/(2jo + 
2j++2). Middle and bottom panels: Occupation numbers nt 
for levels and +. Input parameters: U — 4.0, W = 2.4. For 
the left column, the ratio r = (2j + l)/(2j+ + 1) = 2/3; for 
the right, r — 3.0. 

Results from numerical calculations are plotted in 
Fig. ^ which consists of two columns, each made up of 
three plots. The uppermost panels show the dimension- 
less ratio d+(N) = (e+(N) - e (N))/D. The middle and 



3 



lower panels give, respectively, the occupation numbers 
n and n + . We observe that there are three different 
regimes: in two of them there exist well-defined sp ex- 
citations, and d+ ^ 0, and in the third, the energies of 
the levels and -I- coincide at zero. Passage through the 
three regimes can be regarded as a second-order phase 
transition, with the occupation number n+ treated as an 
order parameter. 

Inserting the above results into Eq. (0}, we find 



E 



M 



E HF (N o =0,N + =N) = -(u-w)(N-N c ) 2 /4 < 



thereby verifying that the M state, having occupation 
numbers < n < 1 for both of the levels and +, has 
lower energy than any HF state. Significantly, the differ- 
ence (|S| is of the same order as a typical shell correction 
5E S in heavy magic nuclei. In such systems, the chemical 
potential fi lies in the large gap between upper filled and 
lower unoccupied sp levels, while in the case of merging 
levels, n is located at the place where the density of states 
attains its maximum. 

The sp levels remain merged until one of them is com- 
pletely filled. If the level fills first, as in the left column 
of Fig. ^ then under further increase of N, quasiparticles 
fill the level +, signaling that the distance d+(N) again 
becomes positive. This behavior resembles the repulsion 
of two levels of the same symmetry in quantum mechan- 
ics, although here one deals with sp levels of different 
symmetry. In the opposite case where level + becomes 
fully occupied before level 0, as in the right column, the 
distance d+(N) becomes negative, and the two levels just 
cross each other at this point. 

In the nuclear many-body problem, both types of sp 
level degeneracy - either initially present or arising in the 
scenario described above - are lifted when pairing corre- 
lations are explicitly involved. The role of -D m i n is played 
by the pairing gap A in the spectrum of sp excitations^ 
To illustrate this situation, we make BCS calculations in 
the above two-level model, under the assurance that re- 
alistic pairing forces are weak enough that the gap value 
remains smaller than the distance between neighboring 
sp levels in magic nuclei. 

This two-level BCS problem is set up and solved as 
follows. First we rewrite the BCS gap equation as 



A = 



gD (Vn (l - n ) + \/n + (l - n+H 



(9) 



In so doing we have followed precedent by introducing 
a common dimensionless pairing matrix element g = 
(2j + 1)Acf /AD, A being a dimensionless pairing con- 
stant. A straightforward derivation, based on the BCS 
identity 4rifc(l - n k ) = A 2 /(e^ + A 2 ), the definition 
e\ = 5E/5n\ with E given by Eq. (QJ, and subtraction of 
one of Eqs. (6) from the other, leads to the key relation 

1 + (U - W)(N - N+) = ^ [R(n+) - R(n )} , (10) 



1), form a closed system determining the occupation 
numbers n Q ,n + and the gap value A. This set of equa- 
tions must be solved numerically; some results are given 
in Fig. [21 The inclusion of pairing correlations does in- 
deed lift the degeneracy of the sp levels. However, the 
value of lowest of the energies Ek of the Bogoliubov quasi- 
particles remains markedly less than D. 

It is instructive to compare the structure of the pairing 
gap A in two cases: when the above shrinkage of the dis- 
tance between the sp levels + and is taken into account, 
and when it is not. In the latter case, A ~ [n\{l — ny)] 1 / 2 
shows two humps with a dip in between^ As seen in 
Fig. |2 the shrinkage effect fills in the dip. This increases 
the part of the ground-state energy associated with pair- 
ing correlations. 
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where R(n k ) = sgn(l - 2n k )y/l/ [4n fc (l - n k )} - 1. and 
(|rU|l . together with the obvious equalities N k = n k (2j k + 
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FIG. 2: Top panels: Pairing gap A (in units of D) plotted ver- 
sus x — N I (2 jo + 2j+ + 2), both accounting for the shrinkage 
of the interlevel distance (solid line) and neglecting it (dashed 
line). Bottom panels: Energies of Bogoliubov quasiparticles 
Ek ~ A/2[rifc(l — nfc)] 1,/2 . Pairing constant: g — 0.3. Other 
input parameters for both columns are the same as in Fig. Q 

Let us now address the case u > w > without pair- 
ing, existing for example in atoms and quantum dots. 
In this case, a pair of particles added to any sp level 
with I always have total angular momentum J =/= 
(Hund's rule), in principle destroying spherical symmetry 
and lifting the m-degeneracy of the sp energies t k m- This 
gives rise to spreading of the levels, the magnitudes of 
which are proportional to u for the level 0, and w for the 
level +. If the interaction function / has long-range char- 
acter, we have u/w ^> 1, and hence the spread of level 
is much larger than that of level +. For the density dp 
associated with the added quasiparticles, we may write 
5p(r) = R 2 nolo (r)J2 m \®li om ( n )\ 2n o m , which is applica- 
ble at least until the crossing of relevant orbitals begins. 
Upon inserting this formula into the relation 5T, = (fSp), 
it is found that the spread does not affect the evolution 
of the centers of gravity e° k = J2m e km/(2jk + 1) of the 
levels, since the isotropic part of 5p has the same form 
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Sop = N R^ olo (r) / An as if the degeneracy of the sp level 
were still in effect. This circumstance is especially impor- 
tant at the stage when the two families of sp levels begin 
to cross each other. Since at u > w > the center of the 
gravity of the level + gets stuck close to the Fermi sur- 
face, our results provide a simple mechanism for pinning 
of the narrow bands in solids to the Fermi surface. 

To exemplify this point, let us consider a model where 
the sp spectrum in local-density approximation (LDA) is 
exhausted by (i) a wide band, which disperses through 
the Fermi surface, and (ii) a narrow one, placed below 
the Fermi surface at a distance D n . We assume that 
only the diagonal matrix element /„„ of the interaction 
function / referring to the narrow band is significant, 
the others being negligible. The shift 6e n in the location 
of the narrow band due to switching on the intraband 
interactions is given by a formula analogous to Eq. (3), 
namely Se n — f nn Pm where p n is the density of the band. 
If the correction Se n exceeds the distance D n then the 
HF scenario calls for the narrow band to be completely 
emptied; but then the shift Se n must vanish. To eliminate 
this inconsistency, only a fraction of the particles leave 
the narrow band, in just the right proportion to equalize 
the chemical potentials of the two bands. The feedback 
mechanism we have described positions the narrow band 
exactly at the Fermi surface, resolving a long-standing 
problem with the LDA scheme. 

In atoms, remnants of an accidental degeneracy of the 
Coulomb problem persist in the formation of electronic 
shells for which the distance between sp levels with dif- 
ferent orbital momenta I is rather small. Recalling that 
matrix elements of the electron-electron interaction are 
quite sensitive to the I value, mergence of definite sp lev- 
els cannot be excluded. To elucidate this situation, one 
needs to analyze the energy functional E = J] efc(0)nfc m + 
\ T. u km,,kim, 1 nk m nk imi , wherein the interaction matrix 
Ukm.ktni! replaces the matrix element J2Jl and summation 
occurs over some states of the last unfilled shell. Results 
from numerical studies of the variational equations gen- 
eralizing Eqs. 0, H = e fe (0) + J2 Wfem.fcimi^feimi, will be 
given elsewhere. 

The new many-body effect uncovered in the forego- 
ing analysis resembles a previously studied phenomenon, 
called fermion condensation, which involves wholesale 
mergence of sp levels in homogeneous Fermi fluids££ In 
any conventional homogeneous Fermi liquid, e.g., liquid 
3 He, the momentum p of an added particle can be as- 



sociated with a certain quasiparticle. Similarly, in most 
spherical odd- A nuclei, the total angular momentum J 
in the ground state is carried by an odd quasiparticle. 
In atomic physics, the electronic configuration of ions of 
elements belonging to the principal groups of the peri- 
odic table repeats that of preceding atoms. From the 
microscopic perspective, in all such "open-shell" systems 
conforming to standard FL theory, the single quasiparti- 
cle term a^^o assumes a special role in the ground-state 
wave function, where represents the ground state of 
a parent system. By contrast, in the case of merging of 
sp levels, the ground-state features a multitude of quasi- 
particle terms and therefore exhibits a more complicated, 
yet more balanced character — as in the comparison of a 
chorus with a dominant soloist. This implication of our 
analysis offers a qualitative explanation of the fact that 
the chemical properties of rare-earth elements differ lit- 
tle, in spite of marked variation in atomic numbers. Such 
an explanation may or may not be at variance with the 
textbook argument^ that the relative squeezing of / and 
d orbitals is responsible for the remarkable similarity. 

In spite of evident commonalities, there is a crucial dif- 
ference between the conditions for the "level-mergence" 
phenomenon in homogeneous Fermi liquids and in finite 
Fermi systems with the degenerate sp levels. In the 
former, the presence of a significant velocity-dependent 
component in the interaction function / is needed to pro- 
mote fermion condensation, while in the latter, sp levels 
can merge even if / is momentum- independent. The rea- 
son for this difference is simple: in the homogeneous case, 
the matrix elements u and w are equal to each other, im- 
plying zero energy gain due to the rearrangement when 
velocity-dependent forces are absent. We point out that 
the study of level mergence in finite systems has the ad- 
vantage of transparency, in that (i) it is free of the com- 
plicated issue of damping sp excitations, and (ii) it gives 
access to the precursor stage of the effect. 

Our exploration of the mergence of single-particle lev- 
els in a finite Fermi system with degeneracy has revealed 
a phenomenon that entails the disappearance of well- 
defined low- lying single-particle excitations, with impor- 
tant implications in diverse physical settings. 
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